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ABSTRACT 


.*v; 


At  a  single  frequency  an  array's  bearing-response  pattern  is  expressed  as  a  weighted 
sum  of  terms  involving  a  steering  matrix  and  the  ergenvalues  and  eigenvectors  of  the 
array's  spectral  density  matrix.  This  point  of  view  is  used  to  gam  an  understanding  of 
the  beampatterns  generated  by  the  Adaptive  Search  and  Track  Array  (ASTA)  processor. 
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AN  EIGENVECTOR  INTERPRETATION  OF  AN 
ARRAY'S  BEARING -RESPONSE  PATTERN 


INTRODUCTION 

Conventional  array  bearing-response  patterns  are  derived  in  terms 
of  the  eigenvector  decomposition  of  the  sampled  acoustic  field's  spectral 
density  matrix.  The  Adaptive  Search  and  Track  Array  (ASTA)  process¬ 
ing  technique  proposed  by  Owsley  1  derives  filters  for  a  multisensor 
array  from  the  components  of  the  eigenvector  corresponding  to  the 
largest  eigenvalue  of  the  spectral  density  matrix.  The  relationship  be¬ 
tween  beam  patterns  resulting  from  this  choice  of  filters  and  bearing- 
response  patterns  due  to  conventional  time -delay  beamforming  provides 
a  basis  for  a  better  understanding  of  both  processing  techniques. 

The  bearing-response  pattern  of  an  n-element  array  of  hydrophones 
can  be  expressed  in  terms  of  the  spectral  density  matrix  of  the  received 
data,  which,  in  turn,  can  be  written  as  a  weighted  sum  of  outer  prod¬ 
ucts.  Each  outer  product  uses  an  eigenvector  of  the  spectral  density 
matrix  and  is  weighted  by  the  corresponding  eigenvalue. 

In  an  acoustic  field  having  equal  power  at  all  hydrophones  and  no 
spatial  coherence,  all  eigenvalues  of  the  spectral  density  matrix  are 
equal.  However,  as  targets  appear,  the  spatial  coherence  of  the  re¬ 
ceived  data  increases,  and  the  eigenvalues,  in  general,  will  no  longer 
be  equal.  Thus,  the  information  contained  in  the  bearing -response  pat¬ 
tern  becomes  concentrated  in  the  terms  of  its  eigenvector  expansion 
corresponding  to  the  larger  eigenvalues.  The  equation  for  a  beam  pat¬ 
tern  of  the  ASTA  processor  is  the  dominant  term  in  the  eigenvector 
expansion  of  the  array's  conventional  bearing-response  pattern. 


DERIVATION  OF  BEARING-RESPONSE  EXPANSION 

Consider  the  conventional  beamformer  of  Fig.  1. 

*N.  L.  Owiley,  ”An  Adaptive  Search  and  Track  Array  (ASTA), "  NUSL  Technical  Memorandum  No. 
2242-166-69,  7  July  1969. 
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Fig.  1  -  Conventional  Beamformet 

Let  the  jxj  (u)j  ,  i  =  1,  2,  .  .  .  ,  n  be  single  frequency,  complex  ran¬ 
dom  processes  at  radian  frequency  u> ,  and  let  the  |e"lwTj|  be  the 
phase  shifts  required  to  steer  the  array  to  some  direction  9  .  The 
|  Hj  (co)|  are  transfer  functions  of  linear  filters  that  may  be  used  to 
"shade"  the  array. 

We  define  the  array's  bearing-response  pattern  as  the  beamformer 
output  power  versus  the  steering  direction  0  .  In  general,  this  pattern 
differs  from  the  array's  farfield  beam  pattern  since  the  latter  is  ob¬ 
tained  by  fixing  the  steering  direction  to  some  angle  9o  and  then  plot¬ 
ting  the  beamformer  output  power  versus  the  angle  of  a  single,  far¬ 
field  source. 

Let  us  define  vectors  x  (w),  y(u>,  9),  and  H(w)  as 


2 


L 


and  a  delay  matrix  D(u,  0)  as 


D  (w,  0) 


Then, 


-jO)T 

e 

— m 


y(co,  6)  =  D(co,  6)  x(oj) 


The  beamformer  output  is  now 

z(ci>,  6)  =  HT(&j)y(«iJ,  6)  «=  y  T(cu,  6)  H(oj)  (2) 

The  bearing-response  function  is  defined  by 

B(oj,  6)  =  E  |  z(cu,  6)  |  2  =  E  [z*  (oj,  0)  z  («,  0)]  , 

where  E  denotes  expected  value.  Substitution  of  Eq.  (2)  into  Eq.  (3) 

yields  ,  *  .i 

B(w,0)  =  e['H  t(w)  y  Kh  0)  yT(w,  0)  H(w)J 

-E  [H^WD'^O)  x^x^d^^HM] 

(4) 

-H*T(w)D*(c4  0)  |  E  [x* (o>)  xT (w)]  |  [  . 

But  the  bracketed  quantity  is  just  the  spectral  density  matrix  of  the  in¬ 
puts;  i.  e.,  Sjj(w)  SlnM 

S  (w)  =  E  [x*(w)  xT(<j)]  =  , 

.  ’  (5) 

S  .(cu)  .  .  .  S  (w) 

n  i  nn  j 

where  Sj.(co)  is  the  cross-spectral  density  function  for  the  ith  and 
jth  inputs.  Thus,  we  can  write 

B(w,  6)  -  H*t(o»)  D*  (oj, 6)  D(ti),  0)H(<i«)  .  (6) 

We  now  take  advantage  of  the  following  expansion  for  S,,(w  )  (See 
the  Appendix): 


S  (w)  =  £  A  (w)  m  (tij)  m* '(oi) 

*  *  _  J  h  n 
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The  |\K(w)!  are  the  eigenvalues  of  S„(oj),  and  the  are 

the  corresponding  normalized  eigenvectors.  Substituting  Eq.  (7)  into 
Eq.  (6)  yields 

B(cu,0)=  1  ^  Ak(w)  [h*t(co)D*(<o,  0)mK(w)]  [mK*T  (w)  D(co,  0)  H(co)]  .  (8) 

But 

•T 

H*T  (cu)  D*  (tu,  0)  mK(cu)  =  [m^,T  (w)  D(to,  d)  H(co)  ]  =  scalar; 

therefore, 

B(cu,  6)  =  1  A  (tij)|  H*tM  D*  (co,  0)  (to)  |  *  .  (9) 

K  =  l  *  * 

Thus,  we  have  expressed  the  array's  single  frequency,  bearing-response 
pattern  as  a  weighted  sum  of  terms,  each  of  which  depends  only  on  a 
single  eigenvector;  the  weights  in  the  sum  are  the  corresponding  eigen¬ 
value  s . 

DISCUSSION  OF  ASTA  BEAM  PATTERNS 
ASTA  has  the  following  structure: 


Tig.  2  -  ASTA  Processor 
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Notice  that  no  steering  delays  are  present.  The  filters  jmK(u>)j  are 
determined  from  the  eigenvector  associated  with  the  largest  eigenvalue 
of  the  array's  spectral  density  matrix;  i.  e.  , 

S  («)  m(ai)  =  A  m(oi)  , 

X*  m  »i 

where  r 

m(o»)  = 

m  (oj) 

n 

It  can  be  shown  2  that  with  these  filters  the  system's  output  power  is 
maximized  under  the  constraint  that 

m*T(ai)  m(ctj)  =  constant 

From  simulation  results,  Owsley3  has  found  that  a  beam  pattern  for 
ASTA  will  have  lobes  in  the  direction  of  the  targets.  Since  the  bearing- 
response  pattern  of  a  conventional  beamformer  also  has  lobes  in  the 
direction  of  the  targets,  one  might  wonder  if  there  is  any  connection 
between  the  two  patterns.  The  answer  is  yes,  as  will  now  be  shown. 

To  compute  a  beam  pattern  for  ASTA,  let  T(w)  be  a  single  fre¬ 
quency,  complex  output  of  a  farfield  source  in  some  direction  0  .  The 
vector  of  received  signals  is  then 

x(cn)  =T(w)D‘l(&i,0)  l=T(w)  D*(u,0)1  , 
where  D(cu,  9)  is  as  defined  above  (Eq.  (1)),  and 


The  beampattern  0(w,0)  is  the  output  power  as  a  function  of  0 
where  0  is  the  direction  of  the  farfield  source;  i.  e.  , 

2 

0)  =  E  |  z(cn,  0)  | 

2J.  N.  Fri.nklin,  Matrix  Theory,  Prentice  Hall,  Inc.,  Lngle  wood  Cliffs,  N,  Y.,  r.'H8,  p.  142. 

3  See  Footnote  1 . 
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But 


Aco,d)  xT  M  m(o>)  =  T ( oj'  m(6J) 


(10) 


As  a  result, 


.  2 

0U..0)  -  Ji:|T(oj)|2[  |  1T  D*  U,0)  mU)|  .  (11) 

Comparing  Eqs.  (9)  and  (11),  we  see  that  for  the  case  of  uniform 
shading  of  the  conventional  processor,  i.  e.  ,  H  (to)  =  H*  (to)  =  1  ,  (3(to,0) 
equals  (to  within  a  scale  factor)  one  term  of  B(lo,  0). 

If  the  acoustic  field  is  of  uniform  power  at  all  hydrophones  and 
possesses  no  spatially  coherent  component,  Sxi(io)  will  be  proportional 
to  the  identity  matrix.  In  this  case  all  eigenvalues  of  Slx(u)  will  be 
equal,  and,  therefore,  each  term  of  (3  (to,  0)  is  weighted  by  the  same 
amount.  A  target,  however,  will  add  a  spatially  coherent  component  to 
the  acoustic  field  and  cause  Sxl(u)  to  deviate  from  a  diagonal  matrix. 
The  eigenvalues  of  Sxx(cu)  will  no  longer  be  equal  in  general;  in  fact, 
certain  terms  of  B(oj,  0)  will  dominate  over  others.  The  ASTA  proc¬ 
essor  provides  us  with  the  term  corresponding  to  \.mgx  ,  i.  e.  ,  the  most 
heavily  weighted  component  of  B(w,  0). 


SUMMAR  * 

It  has  been  shown  that  at  a  single  frequency,  an  array's  bearing-re¬ 
sponse  pattern  can  be  expressed  as  a  weighted  sum  of  terms  involving 
a  steering  matrix  and  an  eigenvector  of  the  array's  spectral  matrix; 
the  weights  in  the  sum  are  the  corresponding  eigenvalues.  This  result 
has  been  used  to  show  that  a  beam  pattern  for  the  ASTA  processor  rep¬ 
resents  the  dominant  component  of  the  conventional  beamformer's 
bearing-response  pattern.  The  results  obtained  here  are  generalizable 
to  the  broadband  case. 
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Appendix 

DERIVATION  OF  EIGENVECTOR  EXPANSION  FOR 
SPECTRAL  DENSITY  MATRIX 

Let  S<(i(u,)  be  a  cr os s - spectral  density  matrix.  Since  S,„(uj)  is 
nonnegative  definite,  it  can  be  diagonalized  by  the  following  unitary 
transformation: 


A ( cu)  -M*T  M (ot)  - 


0 

An  (6j) 


where  the  j  (gj)  J  are  the  eigenvalues  of  S  (u,)  and  M(w)  =  modal 
matrix  of  S  ^  x  ( w)  ;  i.e.,  M(gj)  =  [m  (go)  ,  .  .  .  ,  m  (gj)]  ,  where  the 
(•",  («)[  are  orthonormal  eigenvectors  of  Sx  ^  (gj)  .  "  Solving  for  S  (u>) 
gives 

Cj)  =  M  (o>)  A ( ox)  M  *  T  M 

The  (i,k)  element  of  M(u)  \(u>)  is 

MiKM-  i  L,  (Cl)Y  A  (c>)  5 
,K  q=l  V  q  K  K,< 

where  ( (oj)  ),  =  ith  component  of  mq(u)  and 

j  1  ,  K  =  q 

&Kq  (  0  ,  K  4-  q 

•  T 

The  (i,j)  element  of  M(^)  \(gj)  M  (gj)  is  now 

=  K|,  q|,  (mqM)i  V">  K(l,|)i 

=  |C^,(mKM)i  AK(<  ')  (mK(°',)1  ■ 

Thus,  we  have 

n 

S  (o)  -  1  a.  (gj)  t  -  i.  A  lu)  m.Ui)  m  J.1  (O 

a  *  *  J  |^  _  j  K  K  K 
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At  a  single  frequency  an  array's  bearing -response  pattern  is  expressed  as  a  weighted  sum  of 
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